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Abstract 

A map from the initial conditions to the function and its first spatial derivative evaluated at 
the interface is constructed for the heat equation on finite and infinite domains with n interfaces. 

The existence of this map allows changing the problem at hand from an interface problem to a 
boundary value problem which allows for an alternative to the approach of finding a closed-form 
solution to the interface problem. 


1 Introduction 

Interface problems for partial differential equations (PDEs) are initial boundary value problems for 
which the solution of an equation in one domain prescribes boundary conditions for the equations 
in adjacent domains. In applications, conditions at the interface follow from conservations laws. 
Few interface problems allow for an explicit closed-form solution using classical solution methods. 
Using the Fokas method mm such solutions may be constructed. This has been done in the 
case of the heat equation with n interfaces in infinite, finite, and periodic domains as well as on 
graphs [siiaiisiiiniiii]. The method has also been extended to dispersive problems [13 Hi, and 
higher order problems [3]. These works construct explicit solutions in terms of given initial and 
boundary conditions. The value of the function at the interface is not known. 

In this paper we consider the heat equation with n interfaces on domains of finite and infinite 
extent. The problem of heat conduction in a composite wall is a classical problem in design and 
construction discussed in many excellent texts, see for instance mm- It is usual to restrict to the 
case of walls whose constitutive parts are in perfect thermal contact and have physical properties 
that are constant throughout the material and that are considered to be of infinite extent in the 
directions parallel to the wall. Further, we assume that temperature and heat flux do not vary in 
these directions. In that case, the mathematical model for heat conduction in each wall layer is 
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given by m Chapter 10]; 


, 0 ') _ 


,(i) 


u 


(i) 


(x,t = 0) = 


Xj-l < X <, Xj, 
Xj-l < X <, Xj, 


(la) 

(lb) 


here {x, t) denotes the temperature in the wall layer indexed by (j), aj > 0 is the heat-conduction 
coefficient of the j-th layer, x = xj-i is the left extent of the layer, and x = xj is its right extent. 
The sub-indices denote derivatives with respect to the one-dimensional spatial variable x and the 
temporal variable t. The function Uq\x) is the prescribed initial condition of the system. The 
continuity of the temperature and of its associated heat flux ajUx^ are imposed across the 
interface between layers. In what follows it is convenient to use the quantity aj, defined as the 
positive square root of aj: aj 


la 


j- 


If the layer is either at the far left or far right of the wall, Dirichlet, Neumann, or Robin boundary 
conditions can be imposed on its far left or right boundary respectively, corresponding to prescribing 
“outside” temperature, heat flux, or a combination of these. A derivation of the interface boundary 


conditions is found in m Chapter 1]. It should be noted that the set-up presented in also 
applies to the case of one-dimensional rods in thermal contact. Even for the simple problem of two 
finite walls in thermal contact, the classical approach using separation of variables m can provide 
only an implicit answer. Indeed, the solution obtained in m depends on certain eigenvalues 
defined through a transcendental equation that can be solved only numerically. In contrast, the 
Fokas Method produces an explicit solution formula involving only known quantities. 

The construction of a Dirichlet-to-Neumann map, that is, determining the boundary values that 
are not prescribed in terms of the initial and boundary conditions, is important in the study of 
PDEs and particularly inverse problems [3 EOj • In what follows we construct a similar map between 
the initial values of the PDE and the function (and some number of spatial derivatives) evaluated 
at the interface. This map allows for an alternative to the approach of finding solutions to interface 
problems as presented in earlier papers using the Fokas method by the authors and others. This 
would be most useful in the case where one is interested only in the behavior of solutions at the 
interface. The method presented here can be extended in a straightforward way to many other 
interface problems. To our knowledge, no such maps currently exist. 

Given the initial conditions, one could find the value of the function and its derivatives at 
the interface(s) using these maps. This changes the problem at hand from an interface problem 
to a collection of independent boundary value problems (BVPs). At this point, the BVPs could 
be solved using any number of methods appropriate for the given problem. Each BVP would be 
over-specified, however, by construction it is clear that the corresponding spectral functions are 
admissible |^, i.e. the data is mutually compatible. 


2 The heat equation on an infinite domain with n interfaces 

Consider 

ut = a{x)uxx, ( 2 ) 
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together with the initial condition uo{x) = u{x, 0) and the asymptotic conditions lim| 3 ,|^oc u{x, t) = 
0, where —oo < x < oo, 0 < t < T, and 



X < Xl, 


Xl < X < X2, 

: 

0-n> 

Xn_l < X < Xn, 

2 

X > Xn- 


The restriction lim| 3 ,|^oo u{x,t) = 0 can easily be made more general as in [3]. 

We can rewrite Q as the set of equations 

Xj-i < X < Xj, 0 <t <T, (3) 


for 1 < j < n +1 where xq = —oo and Xn+i = oo. We impose the continuity interface conditions [H 

[ig 


u^^\xj,t) =u^^~^^\xj,t), t > 0, 

a]ui^'>{xj,t) t > 0, 

for 1 < j < n. These interface conditions follow from conservation laws and are fully derived m 
Chapter 1], Since u^^\x,t) is defined on the open interval Xj-i < x < xj, when we write u^^\xj,t) 
we mean lim - u^^\x,t). Similarly, we denote lim + u^^~^^\x,t) by u^^~^^\xj,t). Without loss 

of generality we shift the problem so xi = 0. Using the usual steps of the Fokas method piiniE] 
we have the local relations 


where ujj{k) 


{ajk)'^. These relations are a one-parameter family obtained by rewriting Q. 


t 



T 

1 
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Figure 1: Domains for the application of Green’s Theorem in the case of an infinite domain with n 
interfaces. 


Integrating over the appropriate cells of the domain (see Figure and applying Green’s Theo- 
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rem we find the global relations 


r e-^^^u^o^\x)dx- r 

Jxj — I Jxj — I 

rT 

+ / s) + iku^^\xj, s)) ds 


(5) 


cr 


S) + iku^^\xj-i, s)) ds, 


for 1 < j < n + 1. Define D = {k £ C : Re(wj(A:)) > 0}, Dji = {k £ D : \k\ > R}, and 
= {k £ Dr : Im(fc) > 0} as in Figure 
|x| can become arbitrarily large for j = 1 and j = n + 1, we require k £ C"'' when j = 1 and 
k £ C“ when j = n + 1, in Equation ([^ in order to guarantee that the integrals are defined. For 
2 < J < n, ([^ is valid for k £ C. The dispersion relation ujj{k) = {(Xjk)^ is invariant under the 
symmetry k ^ —k. We supplement the n + 1 global relations above with their evaluation at —k, 
namely, 


2a 


where ii > 0 is an arbitrary finite constant. Since 



Figure 2: (a) The domain for the heat equation, (b) The contour is shown as a red dashed 
line. An application of Cauchy’s Integral Theorem using this contour allows elimination of the 
contribution of terms involving the Fourier transform of the solution. 


PX ' PX' 

0=1 e^'^^ull\x)dx- / ' 

J Xj —I J Xj —I 

pT 

+ / a]e^’^^^+‘^^^’^>iui^\xj,s)-iku^^'>{xj,s))ds ( 6 ) 

Jo 

S) - iku^^\xj-i, s)) ds, 

for 1 < j < n + 1. When j = 1, Q is valid for k £ C“. Similarly, for j = n + 1, Q is valid for 
k £ C"*". For 2 < j < n, Q is valid for all k £ C. Without loss of generality we choose to work 
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with the equations valid in the upper half plane. Define 


9o 


= / e^^u^^\xj,s) ds = / s) ds, 


9i 


(J) 


{uj,t)= e‘^''ui^\xj,s)ds = ^^ e‘^"ui^+^\xj,s)ds, 


Jo 


i^J){k,t)= / ' e-^’^^u^^\x,t)dx, 


u, 


(j) 


sj-i 

Xi 


r^j 

{k) = 

J^.i-1 


e-^’^^u^^\x)dx, 


for 1 < j < n. Using the change of variables k = K/aj on the equation, the global relations 
valid in the upper-half plane are 







( 4,t) -i«> ( =e + 


J 

— inxA_-\ / • 


fj,- 
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-IK 


(7i 


g^^\K\T) + g^^\K\T) 






O-o 


(77 


e'^'Tuin+l) 


— K 


O'n+1 


T -4 


(n+l) 


-/t \ 
<7n+l, 


-]=e^^i^gt\.\T)-^gt\.\T) 

J V(7n+1 ^n+1 

for 2 < j < n. Equation Q can be written as a linear system for the interface values: 

A{k)X{k^,T) = Y{k) + y{K, T), 


where 


A[k , 1) — ygQ , (7g ,..., ^(q , g^ , g^ ,..., g^ 


(n) 


',W 



K 

AY 


'(n) 



. (7n 


— K 
(72 


(«+l) f _ 1 _\\ 

° \crn+l)) 


y{K,T) = r) ,..., uA) r) ,r) ,...,( 


/ —A 


(7b) 


(7c) 


(7d) 


<7rt+l 


and 
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A{k) = 

/ . KX-l 

1 


(7-2 


^2 


—i 
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htx-^ 

^1 



e “^2 


\ 


0'2 


/c 

^2 





<^n 


— IK 

(T2 


e “^2 


iK-i^ 

-g CT„ 

n 



e 


V 


^n-\ ^ 


«^n-l 


—in 

^n 


^ «£ti 

e 


m 

<^n+l 


e 



^^n-l 

<^n 


^ «3;n 
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KXn 
^n + 1 
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The matrix ^(k) consists of four n x n blocks as indicated by the dashed lines. The two blocks 
in the upper half of A{k) are zero except for entries on the main and —1 diagonals. The lower 
two blocks of A{k) only have nonzero entries on the main and +1 diagonals. The matrix A{k) is 
singular for isolated values of k. Asymptotically, for large |k|, the zeros of det(^(K)) lie within a 
strip parallel to the real line m- Since asymptotically there are no zeros in , a sufficiently large 
R may be chosen such that A{k) is nonsingular for every k G and det(^(fi;)) ^ 0. 

Using Cramer’s Rule to solve this system, we have 


ry. _ det(^(^,T)) 

^0 ^ ^ det(^(K)) ’ 

_ det(^+n(At,T)) 
^ ’ det(^(K)) 


(8a) 

(8b) 


where 1 < j < n and Aj{K,T) is the matrix A{k) with the column replaced by T + T- This 
does not give an effective initial-to-interface map because Q depends on the solutions u^^\-,T). 
To eliminate this dependence we multiply Q by Ke~^ * and integrate around D~^, as is typical in 
the construction of Dirichlet-to-Neumann maps [5]. Switching the order of integration we have 




dKds 



dnds 







K det(^j(K, T)) 

det(^(K)) 

Kdet(^j+n(K,T)) 

det(^(K)) 


(9a) 

(9b) 


Using the change of variables il = k? and the classical Fourier transform formula for the delta 
function we have 


Jqd+ det(^(K)) 

uS?>(x,,t) =i / 

' ^^JaD+ det(^(K)) 


(10a) 

(10b) 
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To examine the right-hand-side of (10) we factor the matrix A{k) as A^{k)A^(k) where 

\ 


A^{k) = 


/ —i — xi 

/ e '^1 

1 



1 

1 

1 

1 

— i — Xn 1 

e 1 



7 i—xi 

g ^2 ^ 



1 i — X2 

1 e "^3 

1 
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V 

1 

1 

1 

1 

Xn 

g ‘^n + 1 


is a diagonal matrix. The elements of A^{k) are either 0, 0{k), or decaying exponentially fast for 
K G D^. Hence, 

det(^^(K)) = c{k) = 

for large k in Now, det(^(K)) = c(k) det(^'^(K)) as k —)• oo for k G Similarly, factor 

Aj{K,T) = A^{K)Aj^{K,T)Af{K,T) where Af{K,T) is the 2n x 2n identity matrix with the {j,j) 
component replaced by Then det(^j(K, T)) = det(^^(K)) det(^^(K, T)). Thus, the 

integrand we are considering in is 


/ 

JdDi 


^2^Kdet{Aj{K,T)) 


det(^) 


dK = 


/ 

JdD+ 


^2 Kdet(^j^(K,r)) 


c{k) 


The elements of Aj {k, T) are the same as those in A^ (k) except in the j 
the determinant of A^{k,T) along the column we see that 


th 


dK. 

column. Expanding 


^2^T_t)Hdet{Af{K,T)) _ 
c(k) 


n 


Ci(K) 


-h C£+n(K) e ‘"'^+1 


E 

e=i 


e ’ I —,T] - e 


— K 


<yi+i 


,T]-e 


— K^t — 


IKX£ 


(^+1) 

0 


I — 


— ) 


( 11 ) 


where q = 0{A‘) and C£_|_„ = 0{k) for 1 < .^ < n. The terms involving u^^\-,T), the solutions of 
our equation, are decaying exponentially for k G D^. Thus, by Jordan’s Lemma [I], the integral of 
this term along a closed, bounded curve in C"*" vanishes. In particular we consider the closed curve 
= Cj^+ U where £^+ = dD^ r\{k \ \k\ < C} and Cq = {/c G : \k\ = C}, see Figure 


2b 


Since the integral along Cq vanishes for large C, (0 must vanish since the contour £^+ becomes 
dD~^ as C —)• oo. 

Since the terms involving the elements of y{K, T) evaluate to zero in the solution expression we 
have the solution 









Kdet{Aj{K)) 

det(^(K)) 

Kdet{Aj+niK)) 

detM(K)) 


(12a) 

(12b) 
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for 1 < j < n+1, where Aj{n) is the matrix A{k) with the column replaced by Y[k). Equation[l^ 
is an effective map between the values of the function at the interface and the given initial conditions. 

Remark. Note that since the problem is linear, one could have assumed the initial condition 
was zero for x outside the region X£_i < x < x^. Then, the map would be in terms of just Uq\-). 
Summing over 1 < £ < n + 1 would give the complete map for a general initial condition. 

As an example of a specific initial-to-interface map we consider the equation Q with n = 1. 
Using (12) we have 




Z(T1(72 


7r(ai + (72) JaD+ 


Ke 




( 1 ) 




7r(ui + a2) JdDt 


O’IMq 


2,U1) 


0-1 


^2^?^ ( ^ ) ) dK, 


a^UQ 


— ] + 0-2Uq ^ 
ai 


— 1 1 dAC. 
0-2 


In this case we can deform back to the real line with no pole contributions. Switching the order 
of integration and evaluating the k integral we have 


ud)(o,t) = 


Cri(T2 


2t^Ay/T7{a\ + 02 ) 

1 ( 
y/wt{ai + 0 - 2 ) 



2 2 \ 

-y roo -y \ 

ye^u})\y)dy+ ye^u}, \y) dy ] , (13a) 

Jo 


0 -y^ 


e (y) dy + a 2 I e ^‘""2 (y) dy 1 , 


roo -y ^ ^ 


(13b) 


which is an explicit map from the initial data to the value of the temperature and its associated 
flux at the interface, x = 0. If one allows ai = 02 the problem is simply that of the heat equation 


on the whole line. Equation (13) with ui = 02 is exactly the Green’s Eunction solution of the whole 


line problem evaluated at x = 0 


3 The heat equation on a finite domain with n interfaces 

Consider ([^ on a finite domain, xq < x < x^+i, with the boundary conditions 

/3iu^^\xo,t) + P 2 U^xHxo,t) t>0, (14a) 

/33u("+^^(xn+i,t)+ /34tt^”+^^(x„+i,t) =/2(t), t > 0. (14b) 

As before, we rewrite (© as the set of equations 

=aju'^J^, Xj-i < X < Xj, 0 <t <T, 

for 1 < y < n + 1, subject to the continuity interface conditions 

u^^\xj,t) t > 0, 

fj|4^)(xj,t) =cj|+i4^+^)(xj,t), t > 0, 

for 1 < j < n. Without loss of generality we shift the problem so that xq = 0. 

The following steps are very similar to those presented in the previous section. In what follows 
we give a brief outline of the changes needed to solve on a finite domain. 
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T 

u^^\x,t) 

'u^‘^'>{x, t) 

k 

U^"‘\x, t) 

^(n+l)(^^ t’ 

\ i i 


Xo = 0 Xi X2 ■■ ■ Xn-1 Xn Xn+1 


Figure 3: Domains for the application of Green’s Theorem in the case of a finite domain with n 
interfaces. 


Integrating the local relations @) around the appropriate domain (see Figure]^ and applying 
Green’s Theorem we find the global relations ([^ and their evaluation at —k 0. In contrast to 
Section these 2n + 2 global relations are all valid for k € C. In addition to the definitions in 
Section [ 2 ] we dehne 


JO 

Jo 

Jo 

(w, f) = [ {Xn+i , s) ds, 

Jo 

/m(a;,f)=/' e‘^®/m(s)ds, 


for m = 1,2. Using the change of variables k = K/uj, the n + 1 global relations are 
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Ga 


(15b) 


for 1 < j < n + 1 where we define do = di for convenience. These equations, together with the 
boundary values (14), can be written as a linear system for the interface values 
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where 




= Y^ + y^, 


Y‘^{K,T) = -[-h{K\T),u'i> [-],... Al 


, ^u'o' 1 — ) ^-Mk^^T)] , 

^r). ) \ ^1 J \<^n+l^ 


y^{K, T) = (^0, ,...,«(-+!) ,..., h(”+i) : 0 


T 


and 


A^{k) = 

/ Pi 


<71 


<7l 


<7l 


7^+1 _ 

— m, 


V 


Ps 


P 2 

2 • KXr\ 

(Tr\ —t -^ 

-^e -1 


-CT^ 

• KX-] 

I - - 

e '^1 


—I -!— 

e ^n + l 


_e ‘^rr+l 




The matrix A^ {k) is made up of four (re + 2) x (n + 2) blocks as indicated by the dashed lines. The 
two blocks in the upper half of A^ {k) are zero except for entries on the main and —1 diagonals. 
The lower two blocks of A^{k,) only have entries on the main and +1 diagonals. 

As before we use Cramer’s Rule to solve this system. After multiplying the solutions by Ke~'^ *, 
integrating around D~^, and simplifying as in the previous section we follow a similar process to 
show the terms from y^{ k, T) do not contribute to our solution formula using Jordan’s Lemma and 
Cauchy’s Theorem. One can show that Aj (k, T) can be replaced by AP (k, t) by writing ■ ds 
as Jq • ds + ■ ds and noticing where the function in analytic and decaying. If the boundary 


conditions (14) are time-independent then so is Aj. 


In general, the initial-to-interface map for the heat equation on a finite domain with n interfaces 
is given by 




r ^_^2^ Kdet{Af{K,t)) 

'dD+ iTT det{A^( k)) 


dn, 


JdD+ 77 rdet(A^(fi;)) 

uiHx-1)= f 

^ I ' f7rdet(A^(K)) 


dn, 


dn, 


IdDt 


= [ 

JdDi 


^2iKdet(Af+„+3(K,t)) 


■j+n-|-3v 


iTT det(A^(K)) 


dK. 


(17a) 

(17b) 

(17c) 

(17d) 
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for 1 < j < n + 1, where Aj{ k, i) is the matrix t) with the column replaced by Y^{k, t). 

As an example of a specific initial-to-interface map we consider Q on a finite domain with 
n = 1 with boundary conditions 




t > 0, (18) 

t > 0. (19) 


and zero initial conditions 

0 ) = 0 , 

for 1 < j < n + 1. Using (© we have 


Xj-l < X < Xj, 




10X02 

^ JdD 


ne 


K?, t) sin 


k{xi—X2) 

0-2 


-02f2{K^t)sin{^ 


'r CTlCOS 


k{xi—X2) 

02 


sm 


-^2 cos ( ^ 
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( 1 ) 


(xi,t) =— [ 
TT Ja 




olfx{K?,t) cos 


\ 

n{xi—X2) 

(72 


(7\ 

.2 i (,J1 


sin 

-Cj|/2(K2,t)cOs(^ 


k{xi—X2) 

(72 


( 20 ) 


dK, 


dDo CQg 


k(xi—X2) 


02 


sm 


KXi 

(71 


— (J 2 cos 


KXi 

(71 


sm 


k{xi—X2) 


(72 


dhc. 
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